Distance traveled by random walkers before absorption in a random medium 
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We consider the penetration length I of random walkers diffusing in a medium of perfect or 
imperfect absorbers of number density p. We solve this problem on a lattice and in the continuum 
in all dimensions D, by means of a mean-field renormalization group. For a homogeneous system 
in D > 2, we find that I ~ max K,p~^'^j, where ^ is the absorber density correlation length. The 
cases of D = 1 and D — 2 are also treated. In the presence of long-range correlations, we estimate 
the temporal decay of the density of random walkers not yet absorbed. These results are illustrated 
by exactly solvable toy models, and extensive numerical simulations on directed percolation, where 
the absorbers are the active sites. Finally, we discuss the implications of our results for diffusion 
limited aggregation (DLA), and we propose a more effective method to measure I in DLA clusters. 

I. INTRODUCTION 

The dynamics of random walkers diffusing in the presence of a finite density of perfect absorbers is a rich problem 
which has been widely discussed in the physical and mathematical literature [J. At very large times, the density of 
surviving walkers does not decay exponentially as a simple mean-field argument would predict, but rather behaves as 

n{t) ~ exp -CDp"°+^t"°+^ , (1) 

where p is the absorber density, and Cd is a numerical constant. The physical interpretation is that the process is 
dominated by particles starting in very large absorber-free regions (voids), of linear size L. In D dimensions, these 
regions have a probability of order 

exp(-pL^), (2) 

for small p. In a void of size L, solving the diffusion equation with absorbing conditions on its surface shows that the 
density typically decays as exp(— i/L^). A saddle-point argument then leads to the result of Eq. |^, with the relevant 

regions being of typical size L ^ {t/p)'^+^ , at time t. 

Another important question is the determination of the penetration or screening length I, which measures the 
average distance between the starting point and the absorption point. In the limit of a small density of uniformly 
distributed perfect absorbers of radius a, and for D > 2, a, classic result ^ states that 

However, a simple heuristic argument casts doubts on the validity of Eq. (jS)). In D dimensions, let us consider a 
hypercubic box of linear size L ^ p~^^^ , which is the typical distance between absorbers. We place one absorber of 
radius a in this box, and impose periodic boundary conditions, which is equivalent to copying the box periodically. 
We now release a random walker at time t — 0. To estimate the typical time when the walker will hit the absorber, 
we partition our box in smaller boxes of size a^ . The random walker will be absorbed with a finite probability once it 
has visited most of these small cells, including in particular the one containing the absorber. We define ND{t) as the 
number of different sites visited by a discrete random walker after a time t. Hence, a fair estimate of the absorbing 
time i» is given by 

A^D(i,)~(^) -(P«^)-^ (4) 

During this time, the random walker has traveled a typical distance I given by 

I ^ y/nt^,, (5) 

where k is the diffusion constant. We can use the classical estimates for N^it) [J,y|, which read 

iVi(i) ^ ^, D = l, (6) 



Kf 

Noit) ^ -, D>2. (8) 

Finally, combining Eqs. 14I5I6I7I8|) . we obtain the qualitative estimates 

I ^ p-\ D^l, (9) 



ln(/3a^ 



D = 2, (10) 



l-_D/2 

I ^, D>2. (11) 

In the present work, we will justify on more solid theoretical grounds the results of Eqs. (|9I1UI11|) . In section^ we 
introduce a general formalism in order to compute I and the distribution of distances traveled before absorption. In 
section Hill we define an exact renormalization group for the Green's function on a lattice. This recursion is solved 
using a mean- field (or cavity- like) approximation and our result confirms the estimates of Eqs. (J9I10I1H) . In the limit 
p — > 0, we also address the effect of imperfect absorbers, which becomes relevant in D > 2. In section Hvl we extend 
this approach to the continuum. We compute the distribution of the distances of absorption and its first moments. 
The theory is found to be in excellent agreement with numerical simulations in Z3 = 2 and D = 3. In section we 
apply our renormalization approach to the case of a strongly correlated distribution of absorbers, characterized by a 
power law decay of the absorber density correlation function, c(r) ~ r~", up to the correlation length ^. For a < 2 
and D > 1 (and for a < 1, in D = 1), the penetration length is found to be of the same order as ^, ^ ^ £,. However, 
for a > 2, correlation are weak enough, so that the uncorrelated result of Eqs. 19I1UI11|I is recovered. This result is 
illustrated by the exact solution of the problem in D — 1, and exactly solvable toy models in higher dimensions. In 
addition, we test our ideas on the strongly correlated distribution of active sites in directed percolation simulations 
performed at the critical point, in Z? = 2 and D — 3. Random walkers absorbed by the active sites are found to 
have a screening length I ^ S,, whereas the uncorrelated result is recovered above the critical dimension Dc = 4, when 
correlations become irrelevant. As a conclusion of this section, we extend the result of Eq. ^ to the case of a strongly 
correlated density of absorbers. In section IVII we discuss the determination of the screening length I for diffusion 
limited aggregation (DLA) clusters. We emphasize that the most common method to measure I seems inappropriate 
and we propose an improved scheme. Finally, we give some heuristic arguments inspired by the results of the previous 
sections, leading to the estimate I ^ R, for DLA clusters of gyration radius R. 

II. GENERAL BACKWARD FOKKER-PLANCK FORMALISM 

In a D-dimensional space, we consider a Brownian particle with diffusion constant k. Let V{x) be a positive 
killing field such that if the Brownian particle is a the position x, then in the next time interval dt, it is killed with 
probability V^(x) dt. If the particle is not killed, it just keeps on diffusing. We define P(x, y) as the probability density 
that, starting from x, that the particle's last resting place, i.e. where it is killed or absorbed, is at y. The quantity 
P(x, y) can be calculated by standard backward Fokker-Planck techniques. We consider what happens in the first 
time step dt. One possibility, occurring with probability ^(x) dt, is that the particle is killed where it starts. The 
other possibility is that it is not killed, this with probability 1 — V(x) dt, and so can then make a Brownian jump dB 
having a component in the spatial direction /i, dB^, obeying {dBf^dB^j^i) = 2K5^^'dt. Putting this together gives 

P(x, y) ^ 5(x - y)y (x)dt + (1 - V[^)dt) (P(x + rfB, y)) . (12) 

Expanding to order dt, and taking the expectation value over dB, we obtain 

- «V2p(x, y) + y(x)P(x, y) = <5(x - y)l/(y). (13) 

The solution to Eq. H13|l is given by 

P(x,y)=G(x,y)y(y), (14) 

where G is the Green's function obeying 

- «;V2 G(x, y) + y(x)G(x, y) = <5(x - y). (15) 



Note that the derivation of a probabihty density rather than a probabihty always requires a bit of care and our above 
derivation can be made more rigorous by defining an interval A around the point y then computing the probability 
that the particle is killed in A then taking the limit A — > 0. One can check that P is normalized as follows. Clearly the 
operator acting on G in its defining equation Eq. (|15|l is self-adjoint, which means that G(x, y) = G{y, x). Integrating 
Eq. H15|l over all x, and if the potential V is sufficiently strong, the first term of Eq. (|15|l will give an irrelevant surface 
term. We thus obtain 



/■ 



dxl/(x)G(x,y) = l. (16) 

However G is symmetric which ensures that 

ydyF(y)G(x,y)=ydyP(x,y)-l, (17) 

and which demonstrates the correct normalization of P. To further simplify this problem and reduce everything to 
the study of the Green's function, we write _ff(x, y) = — kVx<5(x — y) and consider V as an operator, V{x,y) = 
(5(x — y)V^(y). Using operator notation, P is given by 

P= [H + VrW ^ [H + Vr^[H + V -H] =1 - [H + Vr^H, (18) 

which reads, in coordinate notation, 

P(x,y)=5(x-y) + AcV2G(x,y). (19) 

Now defining the disorder-averaged values of P and G by p and g respectively, we can write the averaged form of the 
above equation for a translational invariant distribution of absorbers 

p(y) = ^(y) + i^Ksiy)- (20) 

If the disorder is isotropic, we will have p{y) — p{y), where y — |y|, and likewise for g. The disordered averaged 
moments of the distance traveled before the particle is killed are more suitably obtained from p, the Fourier transform 
of p defined as 

p(k) = / dxexp(-ik-x)p(x). (21) 

If we write the Fourier transform of g, g as 

where KE(fc) is given by the one particle irreducible diagrams, we then obtain 

Now, for small k and to leading order, we expect that 

j:{k)^m^ + —k^ + 0{k*), (24) 

K 

where m is the inverse effective screening length of the averaged Green's function and Sk is the renormalization of the 
diffusion constant. The resulting small k expansion for p is 

p{k)^l-^ + 0{k^), (25) 

so if Y is the position where a particle released at the origin is observed, then the disorder averaged second moment 
of Y is simply given by 

2D 

(Y2) = _. (26) 



The characteristic distance from the starting position, /, at which the particle gets absorbed is therefore the same as 
the screening length for the Green's function. We thus define 



l = J{Y^)^l m-\ (27) 

* m 

On a discrete lattice, the same arguments apply and we find that 

p{y) = Sy^o + V^g{y), (28) 

where V^ denotes the lattice Laplacian. The discrete Fourier transform of p is defined as 

p(k)-^exp(-ik-x)p(x), (29) 

X 

and is given by 



p(k) = 1 



2^cos(fc,J -2D 



5(k), (30) 



on a D-dimensional cubic lattice of lattice spacing a = 1. Now, if we write 



^^'')-2i?-2E,cos(M + S(k)' ^^^^ 

we find the lattice result, analogous to Eq. H23fl 

E(k) 
P^^^ ^ 2i^-2E,cos(fc,) + E(k)- (^^) 

The leading order behavior of S must be of the form 

j:{k)Km^ + SKp + 0{k^), (33) 

thus yielding the lattice result 



^2D 



which takes the same form as in the continuum. Interestingly, we find that I is not affected explicitly by the renor- 
malization of the diffusion constant Sk. 

III. MEAN-FIELD RENORMALIZATION GROUP CALCULATION ON A LATTICE 

In this section, we will consider a model on the lattice Z^ (hence the lattice constant is a = 1) with a potential 
VAr(x) given by 

AT AT 

Vn (x) = a ^ (5x,a. =J2u,{x). (35) 

i=l 1=1 

where the a^'s are absorbing sites, with strength A, which arc uniformly and independently (in this first instance) 
distributed among all lattice sites. We denote by V the total number of lattice sites. 

We now estimate the renormalization of the Green's function Gjv for a system with N absorbing sites by the 
addition of an A^ + 1-th absorbing site. A similar method has been introduced in 3] to calculate the effective diffusion 
constant of a tracer particle in a medium composed of randomly placed scatterers. We denote by ajv+i the position 
of the newly added absorber. In operator notation, and introducing the discrete Laplacian H = — kV^, the discrete 
version of Eq. ((T^ for Gat and Gjy+i leads to, 

{H + Vn)Gn = I, (36) 

iH + VN + UN+i)GN+i = /. (37) 



Eliminating H + Vn ~ Gjy i we find 

G^ Gat+i + C/at+iGat+i = J. (38) 

After multiplying by Gat and using the explicit form of Un+i, wc obtain in coordinates notation, 

GAf+iCx, y) = Gjv(x, y) - AGAr(x, aAr+i)GAr+i(aAr+i,y). (39) 

If we set X = a^v+i, we obtain a closed form expression for Gjv+i(ajv+i,y), which can be substituted into Eq. I|39|l . 
This leads to the following exact recurrence: 

r (^ ^r\ ~r (^ .A ^<^Jv(x,ajv+i)GA^(ajv+i,y) 

GA'+i(x,y) - GAr(x,y) , , .^ , x — ■ (,4U) 

Let us first consider the limit where the absorbing sites kill the particle on contact with probability 1. This is 
achieved by taking the limit A -^ +oo and yields 

GAr+i(x,y) ==GA'(x,y) — ^ , (41) 

which describes the exact renormalization of the Green's function due to the addition of a perfect absorber. 

If the perfect absorbers are independently distributed, the position aAr+i is completely uncorrelated with the a^'s, 
for i = 1, ■ ■ ■ , N. Averaging Eq. 141|) over the position aAr+i, but ignoring the correlation between the numerator and 
the denominator in the second term (a mean- field- like approximation), we obtain 

r^ f \ r^ f \ EaG'A'(x,a)GAr(a,y) 

Gat+i X, y = Gat X, y) ^ „ , ^ . (42) 

We now perform the average over the remaining particle positions. Using the statistical translational invariance of 
the system, we find 

^ N / N 1 Ex5w(x;)5w(y-x) 

Taking the Fourier transform of this equation yields 

5iv+i(k)^5iv(k)-i^. (44) 

V.gAf(0) 

Defining the density p — N/V, we obtain a differential equation for the evolution of g{p, k) in this approximation 

op g{p,0) 



where g{p, 0) is self-consistently given by 



r^ dk 
g{p,0)^J ^^^5(p,k). (46) 



The solution to Eq. H45I46|I with the correct boundary conditions is 



AC-l 



^(^'^)-2i^-2E,cos(M + .(p)' ^''^ 



where s{p) obeys 

ds _ r dk 
dp 



_, (27r)^2i^-2E^cos(fc^) + s(p) 

As expected, the defining equation for s{p) does not depend on k for perfect absorbers. 
Now using the fact that s(0) = 0, we can integrate Eq. (|48|l to obtain 



(48) 



'" In 1 + .. .y „„.... Np- m 



{27t)D y 2D-2j:^cos{k^ 

FoT small s, we can use standard results for the integral in Eq. (|48|l |^. The results depend on the dimensionality 
and the lattice structure. 



D = 1: Substituting in the small s behavior of the integral in Eq. H48(l [j|, we find 

ds 



dp 



2^fs. 



For small p, integrating this gives, 



which then leads to 



m^. 


-p^ 


/2 


n/2 


m 


P 



(50) 

(51) 
(52) 



This is clearly the correct scaling in p, which will be recovered when solving exactly the one dimensional case (see 
section HV All . Formally, the result I ^ p^^/^ holds for any dimension \ < D < 2. The physical interpretation 
is clear: the screening length is simply proportional to the mean distance between absorbers, a result which will 
only remain true in the absence of strong correlation between them. 



• D — 2: Here we find 



and integrating this for small p, we obtain 



ds 477 

dp ln(s) ' 

2 ^, 47r/3 



Hence, the screening length is 



m 



HpY 



np 



• D > 2: Here, the leading order is analytic in p |J], and we find 

ds 1 



dp g(0,0)' 



givmg 



I 



l2Dg{0,0) 



The term ^(0,0) depends explicitly on the dimension and the lattice structure, 
dimensional cubic lattice, one has .g(0, 0) sa 0.25... |j|. 



(53) 



(54) 



(55) 



(56) 



(57) 



For example, for the three 



Note that we expect these results to be strongly modified if the absorbers positions are spatially correlated, a 
problem which will be addressed in section IVI 

Let us comment on the disagreement between the present results and the one of Ref. [3 , presented in Eq. ^ . In 
01 , the author first treats the effect of a single absorber on the free Green's function. He then assumes that the total 
correction to S(k) is simply proportional to the number of absorbers. This statement is in fact incorrect, as absorbers 
far away from the introduced random walker should have a negligible contribution. Indeed, the walker should be 
absorbed well before being able to visit the regions where they stand. This is actually the effect of screening, that 
our renormalization approach effectively captures. In addition, Eq. (18) in Ref. 0, which leads to the final result of 
Eq. Q, does not make any sense in the small p limit, and it seems that the opposite non physical limit was in fact 
taken. 

Finally, we can extend our formalism to the case of imperfect absorbers, corresponding to a finite value of A. The 
penetration length I is now expected to depend explicitly on the diffusion constant k. Using Eq. H40|l . we find that 
the Green's function still takes the form of Eq. H47I) . but with s{p) now satisfying 



ds 
dp 



6 



dk 



,{27T)D2D-2j:^cosik^) + s{p) 



(58) 



with 



S= ->0. 
A 



(59) 



In one and two dimensions a random walk is recurrent and visits any position a large number of times. If the density 
of absorbers is small enough, a particle will visit the first absorber encountered a large number of times before making 
an excursion sufficiently far away from this first absorber and visiting a region occupied by a different absorber. This 
means that most particles are absorbed by the first absorber encountered and thus the effect of a finite A should 
become irrelevant in D < 2, when p ^ 0. 
In D = 1, we find the explicit result 



which leads to 



Ss 



2V2S 



VI + 4dp - 1 ■ 



(60) 



(61) 



Hence, for p <C (5 ^, we recover the result of Eq. I|52(l . 
In D = 2, Eq. (jHEl leads to 



ln(s) 

An 



(62) 



/== 



ln(p) ^ 4J 

Tip p ■ 



(63) 



For p <C exp(— 47r5), we recover the result of Eq. l|^ . 

Finally for D > 2, considering imperfect absorbers deeply affects the result of Eq. (|57|l . Indeed, we obtain 



1 = 



l2DigiO,0) + S) 



(64) 



Note that in all dimensions, we find that the penetration length is an increasing function of 5 = j, as physically 
expected. 



IV. THE SCREENING LENGTH IN THE CONTINUUM 



Delta function absorbers in one dimension 



We consider a system where the particle performs continuous Brownian motion in one dimension. It makes perfect 
sense to take an absorbing potential of the form 



N 



V{x) = Xj2S(.x~a,), 



(65) 



corresponding to point-like absorbers. Here, we consider the case where the random walker is absorbed with probability 
one at each absorbing site, that is to say the limit A -^ +00. Without loss of generality, we set k == 1, since the 
final expression of I for perfect absorbers cannot depend on the value of k, whatever the spatial dimension. As in 
the previous section, we apply the cavity approach to calculate the renormalization of the Green's function by the 
addition of an extra absorber at a^+i- We apply again the method of section [llll which leads to the recurrence 
equation 



dp 



5(P,0) ' 



(66) 



where the continuous Fourier transform of g is defined as 

1 f°° 
9{p^x)^--l dk exp{ikx)g{p,k). (67) 

^7^ J-oo 

The solution of Eq. (|66|) is given by 

(68) 



where s satisfies 



^(P)' 



ds ^ , , 

- ^ 2V-S, (69) 



which leads to s = rri^ = p^. Hence, we recover the lattice result 



We also obtain the explicit form 



l^^^^. (70) 

m p 



«') - FT7' '"' 



which gives the distribution of Y, the displacement from the starting position to the point of absorption, to be 

p{y)^^expi-p\y\). (72) 

In fact, in this one dimensional case, the distribution of Y can be computed exactly. A random walker starting 
from X = 0, we denote by b the closest absorbing site to the right and by —a the closest absorbing site to the left. 
Standard results on Brownian motion ^ tell us that the probability of hitting b before —a is 

Pb - ^. (73) 

a + b 

This means that the probability density function for Y before averaging over the disorder is simply 

P(.y) = ^My ~b) + -^6{y + a). (74) 

a^ b a + b 

The disordered averaged density pljj) is now given by 

p{y)^{P{y))a,b. (75) 

where the angled brackets denote the average over the absorber positions a and b. If the absorbers are placed as a 
Poisson point process with rate p, then the probability density function of a and b is Poissonian 

w{x) = pexp(— px). (76) 

This yields 

p{y) = peM-p\y\) r dx ^£^P(Z^ = p T dx (l - ^) exp(-x), (77) 

Jo X + P\y\ Jp\y\ V X J 

showing that in this case the relevant length scale in indeed I ~ p~^. However, we note that the probability density 
function given by the mean-field renormalization method is not exact in one dimension. 

A generalization of this one dimensional model can be constructed as follows. We take the distribution of lengths 
L between the absorbers to be given by q{L). The average interval length is simply related to the density by 



p-' 



/"OO 

= (L) = / dLLq{L). (78) 

"'0 



The probability that the particle starts within an interval of length L has the probability distribution function 

QiL) = ^qW- (79) 

Given that one is in this interval, the position within it is uniformly distributed. This means that we can write 
b = L{1 — U) and a = LU, where U is uniformly distributed over [0, 1]. We thus have 

P{y) = U6{y - L{1 - U)) + (1 - C/)5(y + LU). (80) 

Now performing the disorder average and using the symmetry of the problem, we find 

p{y) ^ / du dLuQ{L)6{\y\-{l-u)L), 
Jo Jo 

\u-^q(-\^). (81) 

Using the explicit form for Q in Eq. (|79|l . we obtain 

/l^ dx (i - M) ,(,) 
Piy) = '" ,oo\ /- . (82) 

Jp dx xq(x) 

As a check, we set q{x) — pexp(~px) and recover the result Eq. 177|l . obtained for the nienioryless distribution of 
Poissonian absorbers. 

The above method of calculation also allows for a straightforward computation of the moments of |y| which can 
be conveniently expressed in terms of the moments of the distribution q{L). We find 

(n + l)(n + 2) (L) 
In particular, we have 

I' = {\y\') - ^^, (84) 

which gives 

l = p-' = {L), (85) 

for uniformly distributed absorbers associated to a Poissonian distribution of intervals. 

B. Absorbing spheres in D > 2 

In this section, we compute the behavior of the distance traveled before absorption in two dimensions and above. 
In this case, we take the absorbers to be spheres of radius a. As before, we restrict ourselves to the limit where the 
density of absorbers p — N/V is small. As in the case of the lattice system, we denote by Gn the Green's function 
in the presence of N absorbers and by Gn+i the Green's function obtained when an extra absorber is placed at the 
point aAT+i which is uniformly distributed in the volume V, independently of the other absorbers. We denote by 
Un+i{^) the potential due to the N + 1-th absorber which takes the form 

f/jv+i(x) = A(5(a- |x-aN+i|), (86) 

where A has the same dimension as n divided by a distance. Upf^i is absorbing with strength A on the surface of the 
sphere of radius a centered at the point a^v+i- In the limit A -~^ +oo, the absorber is a perfect absorber. Combining 
Eq. H15|) written for Gn and Gjv+i, wc can find the equation for Gjv+i in operator notation as 

GjY Gjv+i + f/AT+iGjv+i = /, (87) 
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which can be rewritten as 

Gjv+i(x,y) =G'jv(x,y)- / dzG7v(x,z)L/jv+i(z)Gw+i(z,y)- 

The corresponding equation on the lattice was easy to solve as the potential was given by a delta function. However, 
here we must resort to a further approximation: we assume that we can set GAr(x, z) w GAr(x, ajv+i) and Ga? (z, y) w 
GN{sLN+iTy) in the above integral over the surface of a sphere centered at ajv+i- We obtain 

r< f \ r< f \ Ao-(D)a^"iGAr(x,aAr+i)GAr(aAr+i,y) 

Gn+1 (x, y) = Gat (x, y) . ^ r, jt / \n ( ^ ' (^^) 

1 + J dzC/Ar+i(z)Gjv(z,aAr+i) 

where a{D) is the area of a sphere of unit radius in D dimensions. We cannot make the same approximation for 
the integral in the denominator as the resulting term would be proportional to Gjv(ajv+i, ajv+i) which is finite on 
a lattice, but diverges in the continuous case for D > 2. The finite radius of the absorbing spheres regularizes the 
result. Now, in the limit of large A, we obtain the expression for the renormalization group flow 

r' t \ n f \ o-p)a^"^Gjv(x,ajv+i)GAr(aAr+i,y) 

GAr+i(x, y) = Gjv(x, y) j——- — — ^ , (90) 

where 

/(x, ajv+i) =(5(a- |x-aN+i|). (91) 

Following the same line of arguments as in the lattice case, we find that the disorder averaged Green's function obeys 

dg _ g'^ip.k) 



dp g{a) 



(92) 



where we have used the spherical symmetry of the disorder averaged Green's function. For k = 1 , the solution to this 
equation is 

(93) 



where s obeys 

ds 1 



dp g{a) ' 
and where g{a) has to be computed self-consistently. 
• D ^ 2: In two dimensions, we have 



(94) 



gir) - :^KoiV^3r), (95) 



2tt 

where Ko{x) is the Bessel function of the second kind of order [g. We thus have 

ds 2tt 



dp Kaiy/sa) 



(96) 



Assuming that s is sufficiently small, i.e. for sufficiently small p, we can use the small argument asymptotic 
form of Kq [3, Ko{u) ~ — ln(M), to obtain 

- = ^^ (97) 

dp — In(sa^)' 

where a cut-off of order a^ naturally arises. This equation can be integrated up to the leading order, yielding 




FIG. 1: Plot of (|Y|) (bottom full circles) and I — •\/{Y^) (top full circles) obtained from numerical simulations in D = 2. 
Each point is obtained from the average over at least lO'' random walker trajectories, and we used several samples totalizing 
at least 10* uniformly distributed absorbers. Error bars are much smaller than the size of the circles. We compare these data 
to the functional forms of Eqs. 19911001 (full lines), with the cut-off a^ being the sole fitting parameter. For spheres of radius 
a = 1, we find that both curves are well fitted with an effective value of a^ ~ 2.2. 



We thus find the same functional form as in the lattice case, 
for the lowest order moments: 



After some algebra, we find the following results 



lYD^l,-.!^ 



(99) 



and 



i = ^/(y' 



ln(pa^ 
■Kp 



(100) 



The above calculation actually gives the full distribution for Y. 
for R is given by 



If r = |Y|, then the probability density function 



p{r) = m, rKoimr), 



(101) 



where m is given by Eq. H98() . In the case of imperfect absorbers (finite A), the present results are not affected 
provided that p is small enough {—6/\n{p) <C 1). 

In Fig.^ we plot the two first moments of Y, as found from numerical simulations. We find a perfect agreement 
with Eqs. (|99I100|I . indicating that the theoretically predicted constant prefactors of the p dependence are in 
fact probably exact in the small density limit. Numerically, in order to access to the low density regime, we 
use the following algorithm which has been introduced in the context of DLA ^J . Before performing the next 
move of our Brownian walker, we look for the nearest absorber (by inspection of a grid keeping track of the 
coarse-grained absorber density), say found at the distance d. We then deposit the random walker anywhere on 
the circle of radius d— a centered at its current position. Indeed, the first position where the walker would cross 
the perimeter is uniformly distributed on the circle. The random walker is absorbed when the new distance to 
its nearest absorber is d < a(l -|- e), where e is small enough (typically, we took e < 10~^, and a = 1). Note 
that if we were to actually simulate the random walker trajectory before it reaches the perimeter, the program 
would take a running time larger by a factor d^/di, where dt is the time increment. When accessing to density 
of order p ~ 10~^, and using a time step of order dt ~ lO""^, this factor is of order 10^, which gives an idea of 
the huge gain achieved by using this algorithm. All the simulations performed in this paper use variants of this 
algorithm. 




FIG. 2: Plot (|Y|) (bottom full circles) and I = ^J (Y'^) (top full circles) obtained from numerical simulations in D = 3. We 
compare these data to the functional forms of Eqs. 11U4I105II (straight dotted lines), finding a good agreement at small density. 
Note that a very good fit can be obtained in the entire range of density by using the functional form Ifu = hheory ~ lo, where 
Zo ~ 1, for a = 1 (full lines). 



• _D = 3: In this case, we find 



For small s, Eq. (|M| yields 



?W = 4;^exp(-V^r). 



m — s K 47rap, 



(102) 



(103) 



and we again see that the functional dependence is the same as for the lattice case. The lowest order moments 
are given by 



|Y|) = 



1 

Trpa' 



and 



z = v/rn = y^ 



2'!Tpa 



(104) 



(105) 



(106) 



(107) 



For imperfect absorbers, the result now depends on the diffusion constant 

V 2TTpa V XaJ 

The radial distribution function for Y takes the explicit form 

p{r) — m rexp(— ?nr). 

In Fig. 121 we plot the two first moments of Y, as found from numerical simulations. We find a perfect agreement 
with Eqs. H104I105|I . indicating that these expressions are again probably exact in the small density limit. In 
Fig. 121 we plot the numerical probability density distribution p[r) for p = 10^^, which compares very well with 
our analytical result of Eq. <|107|l . 

• £) > 3: In higher dimensions, we find the same qualitative behavior as in D = 3. In particular, the screening 
length is 



/ = 



2D 



a{D)pa 



D-2 ■ 



(108) 



and is amplified by a factor ^1 + ^ for imperfect absorbers. 




FIG. 3: In three dimensions and for p = 10~ , we plot the numerical probability density distribution of the distance traveled 
before absorption normalized by its average, x = Y/ (Y) (full line), obtained after running a total of 2x10* random walkers. It 
is in very good agreement with the theoretical expression of Eq. Hl()7|l . which leads to p{x) = 4a;exp(— 23;) (dashed line). 



V. THE CASE OF STRONGLY CORRELATED ABSORBERS 
A. Exact result in one dimension 



In this section, we consider the case of a non uniform distribution of absorbers displaying long range density- 
density correlations. In one dimension, such a situation arises naturally when the distribution of intervals q{L) 
between absorbers decays as a power law up to a distance ^, which is much larger than the mean distance between 
absorbers p~^. Hence, we take the typical form 



g(L)~(l-a)-2r^exp(-L/0, 



(109) 



for L larger than a small scale cut-off a (a should not be confused with the size of absorbers, which is irrelevant for 
perfect absorbers in I? = 1; here we consider point-like absorbers). We only consider 



0< a < 1. 



(110) 



The condition a < 1 ensures that q{L) can be normalized even in the limit ^ — > cxd, whereas a > implies that the 
first moment diverges in this limit: 



{L)^p- 



1 



l/^C 



(111) 



Since a < 1, ^ is indeed much larger than the mean distance between absorbers for small p. If the intervals between 
absorbers are drawn independently with the distribution q{L), the density correlation function C{r) — {p{r)p{0)) can 
be computed exactly, through its Laplace transform C{s), which is simply related to the Laplace transform of q{L) 
by the relation |2|] 



C(.) 



m 



1 - qis) 
Hence in real space, the connected density correlation behaves like 



{p{r)p{0))~p' 
Now applying the general result of Eq. (|84|) , we obtain 



■exp(-r/f). 



6(a + 2) 



e 



ai-i/"p-i/". 



(112) 



(113) 



(114) 
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In one dimension, we find that the screening length scales as the correlation length, which is much larger than the 
screening length obtained in the uniform case. In the next subsections, we shall illustrate the fact that this very same 
result should apply in higher dimensions. 

B. A heuristic argument in higher dimensions 

We consider a system where the absorbers are distributed via a physical process which leads to long-range spatial 
correlations between them. Typical examples are given by directed percolation (see next subsection), percolation or 
DLA, the latter being briefly discussed in section IVII In general, the presence of correlations makes the problem 
much more difficult to treat analytically. Here, we present a semi-phenomenological approach based on our mean-field 
renormalization method, which applies in the case where correlations manifest themselves as a clustering phenomena. 
We define the spatial correlation function as 

{p{v)p{Q))^pc{v)+p\ (115) 

where the normalized connected correlation function c(r) is assumed to behave qualitatively as 

^W - '^^^^ (116) 

where £, is the correlation length of the system. We also assume that the system is isotropic. A relationship between 
the density and ^ is found by associating ^ as the characteristic distance where the connected part of the above 
correlation function, the first term in Eq. H115|l . is of the same order as the second term. This gives 

P-C"- (117) 

We implicitly consider the case where ^ is much bigger than the mean distance between absorbers, which scales as 
p-i/D ^ Hence, we will assume from now that 

a<D. (118) 

The regime of small density thus obviously corresponds to a regime where the correlation length is large, for instance 
near a continuous transition. If the correlation is manifested by the formation of clusters, then the typical number of 
absorbers in a cluster is given by 



Mr 



f c{r)dr^^^-", (119) 



which is divergent as ^ — > c», since a < D. If Nc is the number of clusters in the volume V, the total number of 
absorbers is given by A^ ^ Nc£,^~°', which leads to the cluster density 

Pc^pC-^'-r''- (120) 

This last result expresses the fact that the typical distance between clusters of absorbers is of order ^. We thus expect 
to find large empty regions in the system, whose linear size is of order ^, which is, again, much bigger than the mean 
distance between absorbers. In what follows, we will repeat our renormalization calculation but in terms of the cluster 
number. As a starting point, we will use the approximation of Eq. (|89|) . where the potential V is concentrated at the 
center of the clusters and will take the form 

Vir) = X'-^^^, (121) 

which is simply proportional to the mean density of absorbers for a cluster whose center is at r = 0. For large A, the 
flow equation in pc leads to the same functional form for g as before, but the corresponding equation for s — rn? is 
now 

dm? J c{r)r^~^dr 

dpc J c{r)g(r)r^~^dr 

We first consider the case a < 2. In the limit of small m and large ^, we find 

^ ~ ^^-"(m + r ')'-" - ^-'{1 + mO^-". (123) 

Opc 
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where we have used g{r) ^ exp{—mr)/r^^^. Using the relation pc ^ ^~^ , we finally obtain 

K^'-(1 + Kf~"- (124) 

This homogeneous equation admits the obvious solution 

m^C\ (125) 

leading to 

/^^^fli-^/^p-i/"^ (126) 

where we have reintroduced the dependence on the absorber radius a. 

In the case where a > 2, the integral in the denominator of the right-hand side of Eq. (|122|l converges and we 
simply get 

m^ ~ PcC''-'' ^ r" ~ P, (127) 

or 

Z~ai-^/V'/', (128) 

which is the result for uncorrelated absorbers. 

For a = 2 or D = 2, the calculation above leads to logarithmic corrections in the expression of I. Considering the 
crudeness of our argument, we do not believe it is worth detailing the nature of these corrections. 

In conclusion, the present results suggest that the screening length for a correlated system is either the correlation 
length ^ (a < 2) or the screening length found in the uncorrelated case {a > 2). This result can be expressed in a 
synthetic way by 

/~max(c,ai-^/V/'), (129) 

with logarithmic corrections for the uncorrelated result in I? = 2, which were obtained analytically in section HVl 

C. Numerical results for a critical distribution of absorbers arising from directed percolation 

In this section, we will test the ideas presented above by considering a correlated distribution of absorbers generated 
by the active sites remaining at time t, at the critical point of directed percolation j^]- 

Let us briefly introduce directed percolation on a L'-dimensional hypercubic lattice. Lattice sites are empty (inac- 
tive) or occupied by a particle (active). At time t, each site is visited in a parallel dynamics. If the site is occupied, 
the particle is copied on its 2D neighbors with probability p, or removed with probability l—p- If at least one particle 
has been copied on a given site, this site is simply considered as occupied, and empty otherwise. If p is large enough, 
a finite stationary density of particles p establishes itself for large time. On the contrary, for small enough p, the 
density decays exponentially with time. In fact, there exists a critical value pc for which the system is critical: the 
density decreases algebraically, and the spatial (and temporal) correlation length diverges with time. This defines the 
critical exponents 6 and z 

p^r\ i^t^'^ (130) 

and the density correlation function has the typical form introduced in Eq. H115I116|I . with 

a = 5z. (131) 

We have performed extensive numerical simulations of directed percolation at the critical point and considered the 
active sites present at several fixed times. When any of these times is reached, we stop the simulation and launch a 
large number of random walkers which are absorbed by the active sites. We then measure (|Y|) and the screening 
length / = ■\/(Y2), which are found to be proportional. Finally, the directed percolation dynamics is resumed until 
the next sampling time is reached, permitting us to explore systems with smaller and smaller densities, and increasing 
correlation length. We have performed our simulations in Z? = 2 and _D = 3, since our D = 1 result being exact, there 




FIG. 4: Plot of I — ^J (Y'^) obtained from numerical simulations of directed percolation in D = 2 (top dots; 30 samples on 
a 4000^ lattice) and D — ?> (bottom dots; 70 samples on a 300'^ lattice). Each point corresponds to an average over a total 
of 10^ — 10* random walkers, and error bars are smaller than the size of the dots, except for the smallest densities in _D = 2, 
where both are of the same order. The dashed lines have the expected slopes 1/a — {5z)~^ . In both case, the initial curvature 
can be well captured by a two-parameter fit to the functional form I = Ap~ '" + Bp~^' , where the second term corresponds 
to the uncorrelated screening length (full lines). Note that the stronger curvature in D = 2 could as well be due to subleading 
logarithmic corrections as suggested in subsection lV Bl The dotted line represents the slope —1/2, which is the expected result 
in the uncorrelated case, and which holds for D > 4. 



is no doubt that the relation I ~ ^ should be satisfied in this case. For D < 4 {Dc = 4 is the critical dimension above 
which mean-field theory becomes exact), we have a < 2. Interestingly, the equality a — 2 holds exactly in D — 4 and 
above ((5 = 1 and z = 2), so that for D > 4, we recover the result I ~ p~^^^, identical to the uncorrelated case. This 
is comforting, as correlations are known to become irrelevant above the upper critical dimension. For the dimensions 
of interest here, one has 81 



a_D=2 



0.79, 



aD=3 ~ 1-39. 



(132) 



In Fig. ^ we plot I — ^J (Y"^) ~ (|Y|) as a function of p, and find a fair agreement with our prediction I ~ p ^/". 
The numerical data are consistent with a subleading correction of order p~^/^ (the uncorrelated result), although the 
rather strong curvature observed in I? = 2 could be as well ascribed to subleading logarithmic corrections mentioned 
at the end of subsection IV Bl 

We have also measured the average void size A. We pick a point at random in space and determine the radius of 
the largest disk (in Z) = 2) or sphere (in _D = 3) which does not contain any absorber. In _D = 1 (see Eq. (|83ll ). this 
is exactly a measure of ^ since 



A.i<^' 



4 (L) - ¥\ 



l~i. 



(133) 



In higher dimension, we expect this property to hold, since A measures the typical distance between absorber clusters, 
as discussed in subsection IV Bl which was found to be of order ^. In fact, we postulate the more general result 



A ^ max I ^, p 



-l/D 



Fig. |31 illustrates the very convincing linear relation found numerically between I and A, implying I 
have the obvious bound 



(134) 
^. Note that we 



1>X, 



(135) 



as a particle cannot be absorbed in empty regions. 




FIG. 5: Plot of I — •\/(Y^) obtained from numerical simulations of directed percolation in 73 = 2 (bottom dots) and D = 3 
(top dots) as a function of the average void size A (see text), which is expected to be an alternative measure of ^. The fits to 
the linear functional form I = ^A + B are excellent (full lines). The curvature is due to our choice of a log- log plot, in order to 
show that the linear fit works well in the entire range of A. Note that A also seems to share the same statistical fiuctuations as 
^, so that ^ as a function of A in D = 2 is much smoother than when expressed as a function of p (compare with Fig. 2] for 
very small p). 



D. Exact results for some tubular structures in D > 2 



In this section, we consider a specific geometry in D > 2 where the screening length can be computed exactly, 
confirming our general result of Eq. H129|l . 

Let us start by describing the model in D = 2. We consider one dimensional semi-infinite half-lines {x > 0) of 
absorbers separated by a distance L drawn independently with the distribution q{L). Random walkers start from 
a; < 0, and the screening length I is defined as the average depth reached by the walkers before being absorbed. In 
this context, I can be also viewed as a penetration length, similar to the one defined for DLA (see section [^J. A 
random walker penetrating a channel of width L will be absorbed at a depth of order L, the only relevant length 
scale. Actually, the random walker density u in the tube can be computed exactly, by solving the system 



V^M = 0, 
with the absorbing condition for a; > 

u{x, y = 0) = u{x, y = L) = 0, 
and a constant input of walkers at the entrance of the tube 

u{x^O,0<y <L) = 1. 
Using discrete Fourier transform along the direction y, we arrive at the result 

+ 00 



u{x,y) 



-E- 

n=0 



1 



-sm 



1 



(2?! + l)7ry 



exp 



(2n+ l)7ra; 

z 



(136) 
(137) 
(138) 

(139) 



which can also be written as a cumbersome expression involving standard In and arctan functions. The flux of particles 
deposited at the depth x is 



which decays exponentially over the scale L. The average deposition depth is then 

'•+°° Axdx 



m 







Lsmh{^) 



= L. 



(140) 



(141) 




FIG. 6: We consider a semi-infinite tubular system in D = 2 (see text), where the tube width distribution is q{L) 
L e [!,(,]. For a = 3/4, we plot (|Y|) (bottom dots) and I 



L° 



for 



^y(Y'^) (top dots) as a function of the correlation length ^. The 
fits to the linear functional form / = ^4^ + B are excellent (full lines). We also show (L) = p~^ (dashed line), which behaves as 
^", for large ^. Each points is an average over at least 10^ random walkers, involving a total of 10^ absorbing tubes, so that 
the error bars are much smaller than the size of the points. 



If we have an array of such channels, the probability that a walker first enters a channel of width L is 

Q{L) = ^liL). 



(142) 



If we neglect the process where a random walker leaves this first channel to be absorbed in an other one, we find 



/ = 



/ Q{L)l{L)dL = 
Jo 



{L) 



where (L) is the inverse of the absorber density p. If q{L) decreases rapidly, (L^) 

l^p-\ 



and we find 



(143) 



(144) 



In this geometry, the absorbers are highly correlated since they accumulate on lines. Note that the correlation function 
averaged over angle initially decreases as r~^, so that Eq. H144|l is fully consistent with our general result for a — 1. 
Finally, if the width distribution has a power law decay as in Eq. (|109|) , Eq. H143|l leads to 



l^^^p 



-1/q 



(145) 



again in perfect agreement with our general result. 

In Fig. El we present numerical simulations of this two-dimensional tubular system. The data are in perfect 
agreement with Eq. H145|l . showing that the processes involving particles leaving a tube to be absorbed in an other 
one do not affect our general result. 

This model can be generalized in higher dimensions in the following way. We consider linear tubes whose {D — 1)- 
dimensional surface is perfectly absorbing. A (_D — l)-dimensional cut of the system has the structure of a network 
of {D — l)-dimensional "rectangular" cells of edge length Li, L2vj ^i>-ii which are drawn independently using the 
same probability distribution q{L). The absorbers are placed on the {D — 2)-dimensional surface of these cells (see 
Fig. Ejl. The density of absorbers is still 



p={L)-\ 
The solution of the diffusion problem in each tube is simply 



w(a;,2/i,...,yi3_i) 



£1-1 



D-l 



E n 



2n, 

7li ...Tlo — l— J — 1 -^ 



(2nj + l)TTyj 



exp 



'\ 



■^ (2nfc 

k=l 



Ll 
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FIG. 7: In the tubular model introduced in section fVDI the structure consists of linear tubes whose {D — l)-diniensional section 
is a cross product generated by D — 1 linear chains. The absorbers lives on the {D — 2)-diniensional surface of "rectangles" of 
size Li XL2 X ...xLn-i- The figure shows a view from the top of the structure in D = 3, where the sections of the tubes are true 
rectangles. In the model of section IV El the D-dimensional structure is generated by placing the absorbers at the intersection 
of the product of D linear chains. The figure shows the absorbers (dots) in D = 2. 



which decays exponentially on the scale /(L), with 



l(Li) — min L, 

j=l,...,D-l •' 



(148) 



Hence, neglecting again processes where a walker leaves the first tube it enters to be absorbed elsewhere, the average 
penetration length is 



+00 



Jo 



+ OC 



-,_D-2 



dLLq{L) 



Lmin 



If q{L) decays rapidly for L 3> (L), this integral leads to 

However, if q(L) takes the form of Eq. (|109|l . we find 

(•+00 
I - (L)-(^-i) / dLL^L''-^&^^{^Lli)\L''c^^{-Lli)\ 
Jo 



15-2 



(149) 
(150) 

(151) 

(152) 



where we have used the fact that (L) ^ ^" ^ p ^ . We again find Z ^ ^, in agreement with our general argument for 
strongly correlated absorbers. 

E. Exact results for a product distribution of absorbers in D > 2 

In this section, we consider a system for which absorbers are placed at the intersections of the product of D linear 
chains (see Fig.0l. Each chain has its intervals drawn from the distribution q{L). The density of absorbers reads 



P 



(L' 



-D 



(153) 



If the distribution q{L) decreases rapidly on the scale of (L), it is clear that this system will behave similarly to a 
uniform distribution of absorbers, leading to I ^ p^^''^. Hence, we assume that the distribution of intervals is of the 
form 



z(L)^L-2+"/^exp(-L/e). 



(154) 
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We shall see below that a < 2 will lead to / ~ ^, whereas the uncorrelated result is recovered for a > 2. The correlation 
function can be exactly computed, since the density is the cross product of D independent densities 

D 

(p(x)p(0))=[]C(a:,), x=(xi,...,x,3), (155) 

i=l 

where C{x) is the one dimensional correlation function, which can be exactly computed in terms of the Laplace 
transform of q{L) (see Eq. (|112ll '). After performing an average over angle, the correlation is found to qualitatively 
behave as 

(p(r)p(O)) -p"^^ exp(-r/0. (156) 

In the present model, the average void size introduced in section IV CI is defined as the radius of the largest hypercube 
that does not contain any absorbers, averaged over the position of the center of the hypercube. A simple calculation 
leads to the generalization of Eq. (|133|l . If the center x = (xi, ..., xjj) (xi e [0; Li\) is drawn randomly in a cell of size 
Li XL2X ...xix) with all 2^ corners occupied by an absorber, the largest hypercube not containing any absorber has 
a radius 

(i(x) = min mhiixk.Lk — Xk)- (157) 

k=l,...,D 

In order to simplify our calculation, we replace the actual expression of (i(x), by 

(i(x) = min(xfc,Lfc - Xfc), ifc = min Lj, (158) 

j=l,...,D 

which behaves essentially in the same manner as the original expression. Since the average of d(x) over the position 
X of the center is Lfc/4, we finally get 



1 /•+00 r j' + ca 

A=-(i)"'°/ dLrrdn Ll^inliLmin) / dL Lq{L) 

^ Jo VJ L^i„ 



D-1 



(159) 



Interestingly, Eq. (|159|l is the same (up to the 1/4 factor) as the expression obtained for / in the previous subsection 
(see Eq. (jl49|l l except that D — 1 is now replaced by D. Hence, if a < 2, we conclude that 

;>A-f-p-i/". (160) 

In fact, our previous results suggest that ^ ~ A ~ ^ again. Finally, note that in an homogeneous correlated system, 
Eq. H159(l shows that the distribution of the void linear sizes is typically of the form 

QiL)^^eM~L/0, (161) 

which is identical to the void size distribution Q{L) obtained in D = 1 (see Eq. H79|) ). 

Note that in _D > 2, we can obtain a regime where the correlation length is large compared to the average distance 
between absorbers, but is still smaller than the uncorrelated screening length, 

p-^/^«e«p-^/'. (162) 

This regime corresponds to values of a satisfying, 

2<a< D. (163) 

In this case, although there are long-range correlations in the system, our argument of section FV Bl predicts that the 
uncorrelated result I ~ p^^/^ ~ ^"/^ should hold. In Fig.|Hl we measure numerically I as a function of ^ for the model 
studied in the present section, and in D = 3. In perfect agreement with our general result of Eq. (|129|l . we find that 
I ^ ^ for a = 1.5 < 2, whereas I ~ p~^^'^ for a = 2.4 > 2. Note that in both cases, we find A '--^ ^ ^ p^^^^ (not 
shown). 




FIG. 8: Plot of I — ^J (Y'^) as a function of ^ obtained from numerical simulations (full circles) for the model where absorbers 
are placed at the intersections of the cross product of D linear chains in D = 3. For a = 2.4 > 2 (top circles and curves), 
we find an asymptotic slope I ~ ^"'^ ~ P~^ (top dashed line). The fit to the functional form I — A^"''^ + B^ (top full line; 
where the linear correction corresponds to the correlated result) is excellent. As a guide to the eye, we show the line of slope 
unity corresponding to the correlated regime I ~ ,f (top dotted line). For a = 1.5 < 2 (bottom circles and curves), we find that 
I = A^ fits perfectly the data (bottom full line). We also plot a line of slope a/2 = 0.75 (bottom dotted line) to illustrate that 
I ^ ^ grows much faster than the uncorrelated screening length p~^". 



F. Temporal decay of the density of random walkers for a correlated density of absorbers 

In this subsection, we address the generalization of Eq. ^ for a strongly correlated distribution of absorbers in D 
dimensions. We shall adapt to the present problem the usual variational argument [l| leading to Eq. ^. Note that 
the effect of short range clustering on the density of surviving random walkers has been studied quantitatively in J9( . 

Imagine that one releases a random walker of diffusion constant k in an hypercubic region of linear size L not 
containing any absorber (a void). Its survival probability is larger than the survival probability computed assuming 
that the boundary of the void is perfectly absorbing. This probability behaves like P{t) ^ exp(— Kk|i), where k^ is 
the lowest mode of the diffusion equation in the hypercubic domain delimited by the void, with absorbing boundary 
conditions on the surface. All the D coordinates of k^ are equal to ttL~^. By actually solving the diffusion equation 
mentioned above with a uniform initial position of the random walker, we obtain an exact bound of the form 



F(t)>^expf— ^ ^ 



J' 



(164) 



where Ajy is some constant which does not depend on L 

Finally, and after averaging over the void size distribution Q{L), we find an exact bound for the density of surviving 
walkers 



n{t) > Ai 



£""QiL)exp(^-n^^'^dL. 



(165) 



It is commonly conjectured that this kind of bound actually captures the correct asymptotic behavior of n{t) [J. 
Let us now assume a distribution of hypercubic void sizes of the form found in Eq. H161|l . 



Q(L) - ^PiL/O- 



(166) 



In i? = 1, or in the preceding section, we considered a cut-off function F{x) = exp(— x). Here, we will consider the 
more general case 



F{x) — exp(— x^), 



(167) 
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where the value 7 = D is a natural example: L^ j^^ is the volume of the void divided by the correlation volume, an 
equivalent of the term pL^ obtained in the usual uncorrelated case (see Eq. Q). 
Dropping all unimportant numerical constants for the sake of clarity, we obtain 

n(t)-^^°°i-^-^+"exp(-||- (ly) dL. (168) 

Applying a saddle-point argument, and dropping subleading corrections, we finally find that 



n{t) ^ exp I -Ci3 ( ^ ) ) . (169) 



We obtain a stretched exponential decay like in the uncorrelated case, but more importantly, we find that the time 
scale r over which the density decays is now controlled by the correlation length 

r = ^, (170) 

K 

instead of the density p, as obtained in the uncorrelated case 

r=^—. (171) 



VI. PENETRATION LENGTH FOR DLA 

There have been several attempts to measure the screening length of DLA clusters [l3 of gyration radius R 
[llL \12L Ha . I14 . [13 ■ In this section, we show that the methods used so far do not effectively measure I. We will 
propose a theoretical estimate for I as well as a possible numerical method in order to measure it properly. 

First of all, let us briefly mention how DLA clusters are grown in dimension D > 2. One first places a seed at 
the origin, for instance a £)-dimensional sphere of radius a. A random walker (of same radius a) is launched from 
far away and wanders until it touches the seed, whereupon it sticks to it. Then, another random walker is released, 
which sticks upon contact with any of the quenched particles. This process goes on, leading to the formation of fractal 
clusters of dimension Df < D (and presumably, Df > D — 1 10]), defined as 

N^R^f, (172) 

where N is the number of particles in the cluster. We immediately note that for an incoming random walker, the 
points at a distance 2a from the already formed structure acts like absorbers. Hence, it is tempting to introduce 
a screening length I [llL [IJ, measuring how deep the random walkers penetrate the cluster before sticking to it. 
However, contrary to the systems studied so far in this paper, which were infinite and homogeneous, a DLA cluster 
is finite and its density from the center decays in average as 

p(r)-r-(^-^^). (173) 



The determination of I is an important matter, being often an essential ingredient in the theoretical attempts to 
determine Df 

I ^ RiD-Df)/: 

showing that 



ittempts 

determine Df. Although mean- field theories are based on contradicting estimates for I, typically / ~ R^~^f [3ill3 ^^ 
I ^ i?(^^^/)/2 UJ^ they mostly consider that I <C R- However, there is now convincing numerical evidence p^ll3.ll^ 



I - R, (174) 

at least in _D = 2, and probably in D = 3. To understand this result, it is worth mentioning the commonly used 
method to determine I |l3.ll3j| . A cluster of large size R is first grown. Then a large number of independent random 
walkers are released. If one of them touches a particle of the cluster, it simply vanishes, and its last position is 
recorded. The distribution of the positions of absorption is thus obtained and its standard deviation is considered as 
a measure of I (before or after averaging over many clusters having the same number of particles). Although a certain 
consensus seems to exist in the literature, we do not believe that this variance is a faithful measure of I. To see this, 
let us consider the simple example of a perfectly absorbing ellipse in D = 2 of aspect ratio e > 1 and main axis R. Any 
method to measure / should retrieve the trivial result that I = in this case. However, the method presented above 
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would obviously lead to Z ~ i?! It is correct that DLA clusters grown in the continuum are statistically isotropic. 
However, for a given cluster, the fluctuation of the length of the main branches are typically of size R, leading to an 
automatic numerical evaluation oi I ^ R. In other words, even if the random walkers were only absorbed near the 
tip of these branches, physically implying a very small value of I, the current method of estimating I would invariably 
lead to I ^ R. 

As an alternative method (but which might prove as ineffective as the one above), we propose that one should 
measure the distance of penetration from the convex hull of the cluster, and to measure the variance of this length 
before averaging over many clusters. At least, this method gives the correct result for the ellipse (its own convex hull), 
that is, I = 0. 

Finally, inspired by the general conclusions of the previous sections, we would like to give three related arguments 
in favor of the result / ~ R, which we think is correct, although it is argued that the available numerical simulations 
using the method exposed above are not conclusive on this matter. Particles in a DLA cluster are obviously strongly 
correlated, arranging themselves in highly ramified structures. It is thus tempting to apply some of our result to this 
situation: 

• There is no finite correlation length in a DLA cluster except, rather trivially, for the size of the cluster itself. 
For long-range correlations, we have found that I ^^ ^ ^ R. 

• In a DLA cluster, the density correlation function decays as c{r) ^ r^", with a = D — Df. Our study suggests 
that I - p-i/" - R (see Eq. ((T75|) ). 

• Since a DLA cluster is fractal, the largest voids inside it are of typical linear size X ^ R. In homogeneous 
structures, I > A, which suggests again that for DLA, I ~ R. 

VII. CONCLUSION 

In this paper, we have justified analytically the heuristic argument presented in the introduction to estimate the 
screening length in a system of uniformly distributed perfect or imperfect absorbers. Our results in D = 2 and D = 3 
are in excellent agreement with numerical simulations. Even the numerical prefactors of the two first moments are 
surprisingly well described by our approach. As a further check, we found that the theoretical distribution of the 
distances of absorption is in perfect agreement with numerical simulations. 

For correlated absorbers, our analytical approach is not as rigorous as in the uncorrelated case, although we 
supplemented our heuristic argument with the exact solution in D = 1 as well as for two toy models in D > 1, which 
fully confirm our general results. We find that if the density correlation function decays with an exponent a {a < 2 
in D > 2, and a < 1 in D = 1) up to the correlation length ^ ^ p~^^°', then the screening length scales as Z ~ ^ ~ A, 
where A is the average void linear size. These results were confirmed numerically in directed percolation, where the 
active sites play the role of the absorbers. 

Finally, we have argued that for DLA, the penetration or screening length I should be of the same order as the 
linear size i? of a cluster. We have also emphasized that the current numerical method of determining /, although 
finding just I ^ R, could hardly lead to any other result. As a more appropriate method, we propose measuring the 
distance of penetration from the convex hull of the cluster, and to define I as the standard deviation of this length 
before averaging over many clusters. With this definition, the behavior of Z as a function of the cluster radius R is 
certainly a question worth investigating |l8| . 
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